Temperature Dependence of the Dielectric Constant and Resistivity of Diluted 

Magnetic Semiconductors. 
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We study the effect tfiat the ferromagnetic order has on the electrical properties of Diluted Mag- 
netic Semiconductors. We analyze the temperature dependence of the dielectric constant and of 
the resistivity of Gai-ajMuj; As. In our treatment the electronic structure of the semiconductor is 
described by a six band Kohn-Luttinger Hamiltonian, the thermal fluctuations of the Mn magnetic 
moments are treated in the mean field approximation, the carrier-carrier interaction within the ran- 
dom phase approximation, and the transport properties using the relaxation time approximation. 
We find that the Thomas-Fermi length changes near 8% when going from the ferromagnetic to the 
paramagnetic phase. We also find, in good agreement with the experiments, that the resistivity 
changes near 20% when going from zero to the Curie temperature. We explain this change in the 
resistivity in terms of the variation of the Fermi surface and the transport scattering time when 
going from the ferromagnetic phase to the paramagnetic phase. 

PACS numbers: 75.50.Pp, 75.10Lp 
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The incorporation of Mn atoms into III-V semiconduc- 
tors by low-temperature MBE techniques was a big ad- 
vance for the integration of the spin degree of freedom in 
the semiconductor technology. In Gai_a;Mna:As, and for 
Mn concentrations larger than a; '--^ 1%, the magnetic ion 
substitutes a cation, introducing a S* = 5/2 local moment 
and a hole in the valence band of the host semiconduc- 
tor. These compounds present ferromagnetic order, with 
a Curie temperature, Tc near lOOK [l|,l3|j this be- 
havior is generally called carrier induced ferromagnetism 
because the hole carriers mediate the ferromagnetic cou- 
phng between the manganese ions0,|3,IS@ ■ The optimal 
Mn concentration is near x ~ 0.05 and these materials 
are known as Diluted Magnetic Semiconductors (DMS). 
Experimentally!^ it is found that the presence of com- 
pensating defects reduces the number of carriers in the 
system being the density of carriers, p, much smaller than 
that of magnetic ions, c. For the optimal Mn concen- 
tration and not very large carriers densities, most of the 
magnetic properties of the system seems to be reasonable 
well described by using a virtual crystal approximation 
(VGA) Q that neglects thermal and quantum fluctu- 
ations as well as disorder effects. These effects are more 
important as the carrier density in the system becomes 
larger 6, 8, 9] . Post-growth annealing reduces the number 
of compensating defects, increases the density of mobile 
holes in the semiconductor and increases the value of the 
Curie temperature^^ 0|- In the carrier induced ferro- 
magnetism model Tc is proportional to the carrier den- 
sity of states at the Fermi energy, and for this reason the 
annealed samples have a higher Tc than the as-growth 
samples. Another important difference between the an- 
nealed and the as-growth samples is that the later ones 
present a metal-insulator transition at Tc that is absent 
in the former ones0, 1^ the annealed DMS sam- 
ples the resistance increases with temperature (T) near 
25% when going from zero to Tc and remains practically 
constant when T is further increased. As stated above. 



the post-growth annealing process reduces the number 
of defects in the system and the resistivity curves of the 
annealed samples reflect more defect-free intrinsic prop- 
erties of the diluted magnetic semiconductor. 

Most of the theoretical works on DMS study the origin 
of the ferromagnetic order: the effects that the electronic 
properties of the host semiconductor have on the ferro- 
magnetic ground state and the value of Tc- In this paper 
we study the T dependence of the electronic properties 
of the DMS. We analyze how the spin polarization of the 
holes affects the electronic properties of the system and 
from that we obtain the temperature dependence of the 
dielectric constant and of the conductivity of the doped 
semiconductor . 

The main result of our paper is presented in Fig.l. 
There, the T dependence of the electrical resistivity for 
Gai_;j:Mna;As with x ^ 0.05 is plotted for two different 
hole densities. The variation of the resistivity when going 
from zero T to Tc is near 20%, in good agreement with 
experiments in the more intrinsic post growth annealed 
samples ^3 ^3 ■ discuss below, this variation of 

the resistivity with T is due to the dependence of the 
Fermi surface and scattering times on the carrier spin 
polarization and therefore on T . The model we present 
explains the observed behavior of the electrical resistivity 
with temperature. 

In our calculations we describe the electronic structure 
of the DMS using the six band Kohn-Luttinger Hamilto- 
nian and the disorder and thermal effects are treated in a 
mean field approximation. The dielectric constant is cal- 
culated in the Random Phase Approximation (RPA) for- 
malism and the conductivities using the relaxation time 
approximation. 

The system is described by the following Hamiltonian, 



H = Hhoies + J ^ S/ • Si (5(ri - R/) 



(1) 
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Where Hhoies is the part of the Hamiltonian which de- 
scribes the itinerant holes and the last term is the anti- 
ferromagnetic exchange interaction between the spin of 
the Mn^+ ions located at R/ and the spin of the car- 
riers. In order to obtain the temperatm'e dependence of 
the electronic and magnetic properties of the system, we 
minimize the free energy per unit volume, 



T = T 



(2) 



where J-ions is the contribution of the ion spins to the free 
energy and in the mean field description has the form. 
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being h — and ^ the spin polarization of the car- 

riers. Tholes is the free energy of the holes, which is 
obtained in the virtual crystal approximation (VGA) us- 
ing a Luttinger k • p model for describing the carriers. 
In the VGA the average density of states for the real 
system is replaced by that of the average Hamiltonian. 
This approach implies a translational invariant system 
with an effective magnetic field acting on the carrier spins 
Heff = JScm, being m the polarization of the Mn spins. 
In the Luttinger k-p model the wave function of the holes 
in the state (n, k), where n is the subband index and k 
is the wave vector, is expressed as 



(4) 



where |J, mj > are the six T^y valence band wave 
functions. The coefficients a'^^^' and the correspond- 
ing eigenvalues, e„^kj depend on the spin polarization 
of the Mn, and are obtained, from the Luttinger k • p 
Hamiltonian[Tlll3l|. 

By minimizing Eq.Q we obtain the T-dependence 
of the spin polarization of the carriers, ^(T) and of 
the Mn, m{T). In the VGA the Gurie temperature 
has the expression Tc — 2/3cJ^ncr{p) where na-{^J.) is 
the density of states per site and spin at the Fermi 
level [l^ IT3 |. Through the dependence of the Hamilto- 
nian on the spin polarization, we obtain the temperature 
dependence of the eigenvectors, eigenvalues and chemi- 
cal potential. Along this paper we consider always the 
optimal Mn concentration of x=0.05(3j and an exchange 
coupling J=0.060eynm^[l^. Although the properties 
of the system depend on the orientation of the magne- 
tization[? ], this dependence is much smaller than the 
T-dependence and since we are interested in the varia- 
tion of the electronic properties with T we fix the Mn 
spin polarization in the z-direction 

From the eigenvalues and eigenvectors, we calculate 
the dielectric constant that in the RPA has the form, 
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being x(q, ut) the susceptibility. 
X 



(q,^)-2^— — — -T— /(«,k)(i,k+q) ,(6) 



with 
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where np denotes the Fermi-Dirac distribution and ep 
the dielectric constant of the host semiconductor. In the 
GaAs case eo=12.5 

In Fig. 12 we plot the static susceptibility, x(q) = 
x{(l,uj—Q) as a function of the wave vector for 
p—0.02n'm~^, and two different temperatures: T—0 
where the ground state is ferromagnetic with m—1 and 
C=0.777, and T=100K> Tc -57K, where the system is 
paramagnetic, m — ^ —0. The k.p Hamiltonian has cu- 
bic symmetry, and the susceptibility depends not just 
on the absolute value of q, but also on its orientation. 
However, for the range of parameters we are interested 
in, the dependence is almost negligible. We choose al- 
ways the wavector q pointing to the [100] direction. Fur- 
thermore, in the ferromagnetic case, since the magnetic 
anisotropy is small, the dependence of the susceptibility 
on the wave vector orientation is very small and not vis- 
ible in the scale of Fig. |2| The susceptibility depends on 
T mainly at small wave vectors, and this dependence oc- 
curs because the carrier screening ability is proportional 
to the number of states at the Fermi surface, and this is 
smaller in the ferromagnetic phase than in the paramag- 
netic phase. At small wave vectors, only the intraband 
excitations contribute to the susceptibility and the static 
dielectric constant can be written as. 
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and the T dependence of the screening properties can 
be characterized by the Thomas Fermi screening length, 
Xtf = ^/qtf- In Fig. 121 the T dependence of Xtf is 
illustrated. We also show in the inset of Fig. O the T- 
dependence of the ion spin polarization for the same set 
of parameters. The screening length decreases with T 
and changes near 8% when going from the zero T fer- 
romagnetic phase to the paramagnetic phase. Since this 
change is due to the variation of the spin polarization, 
the Thomas- Fermi length remains constant for T higher 
than Tc- 

Finally and in order to compare with experimental in- 
formation, we analyze the transport properties of the sys- 
tem. In the relaxation-time appr oximat ion ' 1 5^ the con- 
ductivity has the form. 
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where r„_k is the elastic scattering time of the electronic 
state (n,k), that has the expression, 
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and Ci is the density of scattering defects with charge 
Qi e. Using equations p 0111(1 for the scattering time we 
have assumed that the main scattering centers are Mn^^ 
acceptors, As-antisite defects, and Mn interstitials, all of 
them have a Coulomb interaction with the carriers [THj. 
At finite temperature the carriers are also scattered by 
short-range spin fLuctuations|l6l IT^. However we have 
neglected the contribution of this spin flip process to the 
elastic scattering time because the Coulomb contribution 
is orders of magnitude larger ^3^^. The conductivity 
of DMS at r=0, for different hole densities and different 
amount of disorder has been calculated recently in Ref. 
[l5l |. Here we study the dependence of the conductivity 
on T. We have evaluated expressions l|^. (|10|) . Hll|l for 
different spin polarization of the holes, which correspond 
to different temperatures. Three factors contribute to 
the variation of the electrical conductivity with T. i) 
A change of the thermal Fermi distribution of the car- 
riers. In the range of temperature of interest this ef- 
fect produces only a small variation of the electronic 
properties 19j. ii) The change of T produces a change 
in the dielectric function Fig|21 and therefore a change in 
the matrix element Ec ljlll) . However this effect is very 
weak because in the transport scattering time the for- 
ward scattering is suppressed by the (1 — cos^k^k') term 
and the more important matrix elements have associated 
a large transfer of momentum for which the electrical 
susceptibility is practically T independent, see Fig|21 iii) 
The change from the ferromagnetic to the paramagnetic 
ground state produces a reduction of the majority spin 
Fermi surface area. For that reason the wave vector dif- 
ference that appear in the matrix element Eq. (|ll|l is big- 
ger in the ferromagnetic phase than in the paramagnetic 
phase, and therefore, the scattering time and the con- 
ductivity are larger in the ferromagnetic phase than in 
the paramagnetic phase. This is the main contribution 
to the variation of the conductivity with T. 

In Fig.lQJ we plot the resistivity of Gai_2;Mn2;As, 
for x=0.05, J—60'meVnnT' and two different densities, 
p=0.2nm-3 (Tc = 57K) and p^O.Gnm-^ (Tc = 125K). 
We normalize the resistivity to the value of the resistivity 
in the paramagnetic phase. As the temperatures studied 
are always much smaller than the Fermi temperature, 
in the paramagnetic phase the resistivity is practically 
T independent. As reported in ref.([T3) in the ferro- 
magnetic phase there is an anisotropic effect, and the 
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FIG. 1: Resistivity, normalized to its value in the paramag- 
netic phase, as a function of T for Gai-^Mnj; As, with a;=0.05, 
J—Q.Q&eVnm? and two different hole densities, p=Q.2nm~'^ 
and p=Q.Qnni~'^ . 



resistivity is different for directions parallel and perpen- 
dicular to the magnetization. The anisotropy is 2.5 % 
for the p=0.2nm~^ case and smaller than 0.5% for the 
p—O.Gnm'^ case. More interesting, for both hole densi- 
ties there is a bing change, near 20%, in the value of the 
resistivity when going from zero T to Tc- This change is 
very similar to the experimentally observed (lOl ITI| in the 
post growth annealed samples. Also the overall shape of 
the resistivity versus T curve is similar to the experimen- 
tal one. 

Therefore we conclude that our model, that is based in 
the use of a)the six band k-p model for the band structure 
of the DMS, b)the mean field approximation for describ- 
ing the thermal fluctuations, c)the RPA for the dielectric 
constant and d)the relaxation-time approximation for the 
resistivity, describe appropriately the electric properties 
of DMS and their dependence on T. Our theory is able 
to describe adequa tely the experimental dependence of 
the resistivity on T p^lTH . 

From our calculations we know that the variation of 
the resistivity when going form the ferromagnetic phase 
to the paramagnetic phase is due to the different Fermi 
surfaces that have these phases. The wave vector trans- 
ferred in a scattering event across the Fermi surface is 
larger in the ferromagnetic phase than in the paramag- 
netic phase and this produces that the resistivity in the 
paramagnetic phase is larger than in the ferromagnetic 
phase. 
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FIG. 2: Static electrical susceptibility as a function of the 
wave vector for a Mn concentration a;=0.05, a hole den- 
sity p=0.2nm~^ , an exchange coupling J=0.06eVnm^ and 
two temperatures: T—0 (Ferromagnetic ground state) and 
r=100K (Paramagnetic ground state), a is the FCC lattice 
parameter of the host semiconductor. 
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FIG. 3: Temperature dependence of the Thomas-Fermi 
length for a Mn concentration a;=0.05, a hole density 
p=0.2nm~^ , and exchange coupling J=0.06eVnm^. The in- 
set illustrates the T dependence of the Mn spin polarization. 
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